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Abstract 

Recent studies of several multi-body meson decays have revealed that the final 
states are dominantly CP-even. However, the small value of the width difference 
between the two physical eigenstates of the system indicates that the total 

widths of decays to CP-even and CP-odd final states should be the same to within 
about a percent. The known contributions to the width difference from hadronic 
decays are discussed, and it is shown that an apparent excess of quasi-CP- 
even modes is balanced, within current uncertainty, by interference effects in quasi¬ 
flavour-specific decays. Decay modes which may significantly affect the picture with 
improved measurements are considered. 


ii 


1 Introduction 


Oscillations between the flavour eigenstates of the neutral D meson result in physical 
states with distinct masses and widths. Parameters that quantify the mass and width 
differences are given by 


xd = {m 2 -mi)/TD and yo = {T 2 -Ti)/{2Td) , ( 1 ) 

where mi( 2 ) and ri( 2 ) are the masses and widths of the physical eigenstates and Pd = 
(Pi+P 2)/2 is the average width. These are now known to good precision from studies of the 
evolution with decay time of the rates for certain D decays, the latest world averages [1] 
giving 

XD = (4.1 tli) X 10-3 , yo = (6.31°^) x lO'^ . (2) 

The branching fractions for decays to CP-even and CP-odd final states are related 
to the reduced width difference j/d, in the limit of CP conservation, by 

r(pO ^ CP-even) - r(D^ CP-odd) B(D^ CP-even) - B(D° CP-odd) , , 

® =-2f^-=-2-■ 

In Eq. ([3]) the widths and branching fractions to CP-even and CP-odd states include not 
only decays to CP eigenstates, but any decays with net CP content [2]. This can be 
quantified by the fractional CP-even content P+ (P+ = 1 corresponds to pure CP-even, 
P+ = 0 to pure CP-odd). Decays with net CP content include quasi-CP eigenstates (he. 
self-conjugate final states with P+ 7 ^ 1) and also quasi-flavour-speciflc decays. Considering 
these different types of decays, the relation for y^ can conveniently be expressed 

VD = \Y.{2F^i-mD° ^ h), (4) 

i 

where the sum is over all distinct decays to hadronic final states /j, each with net CP 
content 2 P+j — 1 . Semileptonic and other nonhadronic decays do not contribute as they 
are either assumed to be flavour-specific, hence 2P+ — 1 = 0, or have branching fractions 
that are negligibly small. 

Measurements of the P+ values of various decay modes have recently become available, 
from studies of samples of quantum-correlated -0(3770) ^ DD decays. These include 
direct measurements as well as determinations of quantities that can be translated to 
provide information on P+. It is therefore of interest, and is the main purpose of this 
paper, to evaluate Eq. (jl]) from experimental data, and to compare the value obtained 
with that from Eq. ([2|). Throughout the paper CP conservation in mixing and 

decay, which is known to be a good approximation, is assumed. 

In addition to the intrinsic value of such a survey, the outcome may be useful to identify 
additional D decay modes which could be used to study charm mixing and CP violation, 
and also for determinations of the angle 7 of the Unitarity Triangle formed from elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [3111]. Decay modes 
with high net CP content are particularly useful to probe these phenomena. Many of 
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the CP eigenstate modes with largest branching fractions, however, give decay topologies 
{e.g. KgTi^) withont any charged track originating from the D decay vertex, and therefore 
cannot be nsed to stndy decay-time-dependent effects. If decay modes with high net CP 
content and experimentally accessible topologies can be identihed, they conld be nsed to 
improve the precision of charm mixing and CP violation parameters [5]. 

The CKM angle 7 can be determined from B —)■ DK decays with negligible theoretical 
nncertainty when the nentral D meson is reconstrncted in hnal states to which both 
and can decay. Methods have been proposed, and analyses performed, with D decays 
to CP eigenstates (referred to as the GLW method) [6l[7], donbly-Cabibbo-snppressed 
decays (ADS) [Hlin] and self-conjngate mnltibody decays (GGSZ) [TUl - lT^ . Ultimately, 
the best precision on 7 is obtained by combining resnlts from all of these methods and 
others [Tdl - fTH] . and it is important to continne to improve the precision in all of them. 

In case the hnal state of a mnltibody D decay is dominated by a particnlar CP- 
eigenvalne, it can be nsed in a “qnasi-GLW” analysis HZ]. First qnasi-GLW analyses 
of B ^ DK with D 7 r+ 7 r“ 7 r° and D have recently been reported by 

LHCb [IB]. If additional decay modes with high net CP content can be identihed, they 
conld also be nsed to improve the precision on 7. Similarly, mnltibody qnasi-havonr- 
specihc decays can be nsed in a “qnasi-ADS” analysis [91119], as has been done with 
D —)■ [1^1201121] and [22] decays; in this case the sensitivity depends 

on a “coherence factor” [19] that is also related to the net CP content of the hnal state, as 
will be shown. Therefore modes with high net CP content are also very nsefnl to improve 
the sensitivity to 7. 

The remainder of the paper is organised as follows. In Sec. O the cnrrent knowledge 
of charm decays to CP eigenstates and (qnasi-)CP eigenstates is reviewed, while Sec. [3] 
contains a similar discnssion of qnasi-havonr-specihc channels. In Sec. HI channels poten¬ 
tially worthy of fntnre experimental investigation are identihed. A brief snmmary is given 
in Sec. El 

2 Charm decays to (quasi-)CP eigenstates 

The cnrrent knowledge [23] of branching fractions of hadronic decay modes to CP 
eigenstates or qnasi-CP eigenstates is given in Table [T] The corresponding valnes of 
2A+ — 1 are also given. Snmming all these valnes according to Eq. (Hj) gives 

^n-qcp ^ Q ^ ^ 2) X 10"^ , (5) 

where the notation indicates that this valne is determined from the experimental 

knowledge of the branching fractions of decays to qnasi-CP eigenstates only. It is not a 
trne estimate of yn as there are additional modes, inclnding those discnssed in Secs. |3]and|l| 
that are not incinded in the snm. Note that the nncertainty on y^^^^ is dominated by the 
least precisely measnred modes, with the largest contribntions coming from nncertainties 
on the branching fractions of D decays to and hnal states. 

Table [H inclndes the mnltibody D decay modes for which explicit determinations of 
P+ have been made [T71[2^ : 
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• ^ 7r+7r-7r°: F+ = 0.973 ± 0.017, 

• ^ K+K-tt^: F+ = 0.732 ± 0.055 , 

• ^ 7r+7r-7r+7r-: F+ = 0.737 ± 0.028 . 

The dominance of CP-even in —>■ decays is particularly striking. It should be 

noted that this feature has been previously discussed in terms of isospin [25] and of the 
impact on the determination of 7 from Dalitz plot analysis of P —> DK, D —)■ 
decays [26]. A potential explanation in the context of flavour-SU(3) and factorisation has 
been discussed m but there is as-yet no fundamental understanding of the CP-even 
dominance in these decays. 

For some other multibody D decay modes, although there is no direct determination 
of P+, it is possible to obtain constraints from published information. For —> Kg7i'~^7r~ 
or KgK^K~ decays, this can be done in one of two ways: (i) from the ratios of double- 
tagged and single-tagged yields for each CP-eigenstate; (ii) from the relation between P+ 
and the factors Ki and Q ini, 

n = ^ 1 (Ki + K-i + 2c, /yTrl") . (6) 

i>0 

where Ki is the fraction of flavour-tagged D decays that fall into Dalitz plot bin i 
— 1 ); D is decay-rate weighted average of the cosine of the strong 
phase difference between P° and P° decays in bin i [Tnll32l l33]. and bins i and —i are 
related by symmetry under charge conjugation of the hnal state. The values quoted in 
Table [T] are determined with the second method, which should be more precise. The same 
approach can also be used to determine for —)■ K^n~^7r~, 2P+ — 1 = —0.288 ± 0.036 

and for P° —> K^K^K~, 2P+ — 1 = —0.136 ± 0.088; since the branching fractions for 
these modes have not been measured they are not included in Table [T] (but see discussion 
below). For the mode, only approach (i) is possible; an estimate of 

P+ could in principle be made from the amplitude model for that decay [31], but the 
evaluation of uncertainties would be difficult. 

The majority of CP-odd hnal states are of the form {K = 7 r°, 7 , p', cu). For 

each of these there is a corresponding CP-even hnal state, K^h^. Among the decays 
with mesons in the hnal state, only the branching fraction for has been 

measured [3l]. The result dihers from that for P° —> by about 4tan^6'c ~ 20 % 

of its absolute value {9c is the Cabibbo angle) due to interference between the Cabibbo- 
favoured and doubly-Cabibbo-suppressed amplitudes [55]. Similar relations are expected 
between the other K^Ji^ branching fractions where K is a pseudoscalar {'r],r]') [361155] . 
though somewhat diherent ehects are possible in other cases [39]. Comparisons of the 
2P+ — 1 values for JFg and K^^K^K~ discussed above also indicate approximately 

opposite CP content between the hnal states containing iCg and mesons. 

Thus, for many of the CP-odd hnal states, there are expected to be corresponding 
CP-even hnal states with branching fractions of comparable magnitude. The converse is 
not true: for the majority of CP-even hnal states (apart from the states), there 
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Table 1: Branching fractions of hadronic decay modes to CP eigenstates or quasi-CP eigen¬ 
states. Values are taken from Ref. [23] unless otherwise specified. The particles vr, K, rj, r( and a; 
are considered as stable; decays involving other particles, such as p or </> mesons, are accounted 
for under the corresponding multibody final state. Effects of CP violation in the neutral kaon 
system are at the level of 0(10“^) and are negligible. 


Decay mode 

2P+ - 1 

P/10-3 


K+K- 

CP-even 

1 

3.96 

±0.08 


KIKI 

1 

0.17 

±0.04 


povro 

1 

10.0 

±0.7 


7r’''7r 

1 

1.402 

±0.026 


7r°7r° 

1 

0.820 

±0.035 


7r°?7 

1 

0.68 

±0.06 

[231128] 


1 

0.90 

±0.14 


TT^U 

1 

0.11 

±0.04 


?7?7 

1 

1.67 

± 0.20 


7]7]' 

1 

1.05 

±0.26 


KIKIKI 

1 

0.91 

±0.13 


KgTT^TT^ 

1 

9.1 

± 1.1 


KgrjTT^ 

1 

5.5 

± 1.1 


7r+7r“7r° 

Mostly CP-even 
0.946 ±0.034 [21] 

14.3 

± 0.6 


K+K-tt^ 

0.464 ±0.110 [21] 

3.29 

±0.14 


-1- 

1 

-1- 

1 

0.474 ± 0.056 [21] 

7.42 

± 0.21 


/CgTr+TT 

Approximately CP-neutral 
0.112 ±0.024 rann] 28.3 

± 2.0 


K^K+K- 

0.194 ±0.064 [29l[30] 

4.60 

±0.16 

[ 231128 ] 

K+K-Tl+Tl- 

0.14 ±0.21 [H] 

2.43 

± 0.12 


Pg°7r0 

CP-odd 

-1 

11.9 

±0.4 


Kin 

-1 

4.79 

±0.30 


K'W 

-1 

9.4 

±0.5 


A> 

-1 

11.1 

± 0.6 


o 

o 

o 

-1 

< 

0.4 



-1 

< 0.6 



is no corresponding CP-odd final state. Essentially, this is due to the fact that singly- 
Cabibbo-suppressed transitions appear to result predominantly in CP-even hnal states, 
while Cabibbo-favoured transitions can give both CP-even and CP-odd states. Excep¬ 
tions arise from the three-body CP-even states Pg 7 r° 7 r°, KgTr^r] and K^K^Kg, for which 
corresponding CP-odd states K^n^r] and K^K^Kg exist. (Although numerically 

insignihcant, it is interesting to note that there is no CP-odd counterpart to the 
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mode, as the symmetry of the wavefunction prevents —)■ decays.) 

This raises a subtle issue in the evaluation of Eq. (|4]). While it is valid to sum all modes 
for which experimental results are available, this could be argued to be biasing, since more 
data are available for modes containing than mesons. Considered together, the net 
CP content can be expected to cancel to some extent, although the residual contributions 
may still be important. For comparison, if all such modes containing mesons in the 
hnal state are excluded from the calculation of a larger value of (14.7 ± 0.5) x 10“^ 

is obtained. 


3 Quasi-flavour-specific hadronic charm decays 

Decay modes such as K^tt^ have contributions from both Cabibbo-favoured and 

doubly-Cabibbo-suppressed amplitudes. The interference between these two amplitudes 
leads to an asymmetry in the widths for CP-even and CP-odd decays [401 - 142] . 

aCP _ ^ {Dcp- f) — ^ {Dcp+ f) _ Up — Rf^f cos Sf 

f - + 1 + Pws/ ■ 

Here the hnal state / includes both conjugate states (e.g. f = Vf is the average 

ratio of magnitudes of, and —Sf is the average strong phase between, the suppressed and 
favoured amplitudes, P/ is the coherence factor, 0 < Pj < 1, that quantihes the dilution 
due to integrating over the phase space (for a two-body decay is the suppressed- 

to-favoured amplitude ratio and P/ = 1), and Pws / is the decay time integrated ratio of 
wrong-sign to right-sign rates, including effects from both the doubly-Cabibbo-suppressed 
amplitude and charm mixing. Precise dehnitions of these quantities can be found, for 
example, in Ref. [IH]. Note that alternative conventions for the dehnition of Sf can be 
found in the literature; in particular, Sf ^ is also widely used, e.g. in Ref. [43]. Since 
a positive sign for in Eq. ([7]) indicates a larger width for CP-odd decays, neglecting 
small corrections from yp and Pws, leads to the relation 

2F+f-l = -A^^ = 2RfrfCosSf. (8) 

The asymmetry for —)■ K^tt^ decays has been measured to be = (12.7 ± 

1.3 ±0.7) % [13] (see also Refs. [44] - l45] b Similar quantities for other quasi-flavour-specihc 
decay modes can also be obtained from published results. Measurements of the relevant 
properties of P° —)■ and decays have been performed |47l|48], where 

the reported quantity Acp is equivalent to A'^^. These values are given in Table [2j It 
should be noted that the signihcantly non-zero value of Acp for P° ^ P^tt^tt+tt” decays 
is in tension with the small value of the coherence parameter for these decays also reported 
in Ref. [IS] . This arises as the latter uses information from all different tags reconstructed 
in the '0(3770) DD events, while Acp is evaluated using CP tags only. This shows 
that a more precise determination of 2P+ — 1 for this mode is possible using Eq. o or (]8]) 
and all available information on the relevant parameters. Improved measurements with 
larger -0(3770) data sets are clearly well motivated. 
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Table 2: Branching fractions of hadronic decay modes to quasi-flavour-specific final states. 
Values are taken from Ref. [23] unless otherwise specified. 


Decay mode 

2P+ - 1 


P/10-=^ 

Cabibbo-favoured / doubly-Cabibbo-suppressed 


-0.127 ±0.015 

m 

38.8 ±0.5 


-0.084 ±0.028 

m 

139 ±5 

K^7r^7r^7r~ 

-0.119 ±0.027 


80.8 ±0.20 


Singly-Cabibbo-suppressed 
0.65 ±0.14 [19] 

5.6 ±0.6 


For completeness, it should be noted that there are also singly-Cabibbo-suppressed 
decays to non-self-conjugate final states, an example of which is This 

decay has been studied |19] with results that indicate significant net CP-even content. 
The relation — k~) / {k'^ + k~), which is valid neglecting charm mixing, can 

be used to together with the values of reported in Ref. [19] to obtain = 

—0.65 ± 0.14 (small possible correlations have been neglected). Alternatively, the right- 
hand side of Eq. (|7|) can be evaluated using quantities determined using all tag information. 
(In this case it is not valid to assume rj 1, as done in Eq. ([ 8 |), though charm mixing 
effects can still be neglected.) This approach gives a consistent value of A^^= 
—0.69±0.08. It would be interesting to see if the dominance of CP-even in —)■ Pg 

decays is confirmed with studies of larger '0(3770) —> DD decays, or if it is consistent 
with the recently-obtained amplitude models for these decays [50] • (To obtain an estimate 
for 6 from amplitude models of P° —?■ JFgiF^ 7 r“ decays requires knowledge of the relative 
phase between the amplitude in the two distinct final states. This quantity, specifically 
the phase between the amplitudes for P° ^ K*^K~ and K*~ decays, can be and has 
been determined from Dalitz plot analysis of decays [5THS3]-) 

Summing the values from Table [2] according to Eq. (|1|) gives = (—11.3 ± 2.3) x 

10 “^, where the notation indicates that this determination comes from the mea¬ 

sured quasi-flavour-specific modes. If modes with final states containing mesons are 
excluded, for the reasons discussed in Sec. O this becomes (—13.1 ±2.3) x 10“^. The 
uncertainty is dominated by the imprecise knowledge of 2 P+ — 1 for P° —>• and 

A^ 7 r^ 7 r’'' 7 r“ decays. It is striking that all of the D —> and P^Tr^Tr+Tr” 

decays are net CP-odd. 

4 Possible additional modes with net CP content 

Using Eq. (|1|) and the resulting relation yj^ = allows a 

prediction of the net CP of the decay modes that have not been accounted for, due to 
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absence of experimental information, 


missing 

Vd 


i ( 2 F+,. - 1 )B(B" ^/i) 

i (missing) 


(6.6 ±2.6) X 10■^ 


(9) 


or (4.7 ± 2.3) x 10“^ if final states containing mesons are exclnded. Since, the Particle 

Data Gronp reports that (38.2 ± 1.4) % of decays are as-yet nnacconnted for, this 

indicates that these missing modes either have low valnes of 2F+ — 1, or that there is a 
balance between net CP-even and CP-odd modes leading to the low valne of Eq. dH]). (For 
comparison, the modes included in Tables [T] and [2] account for about 40% of the total 
width - the remainder are predominantly from semileptonic decays.) 


Table 3: Hadronic decay modes that could be of interest with regard to their CP content. The 
branching fractions are given where known [23]. The branching fraction for —>• 

{K^7r^Tr~^7r~Tr^) includes contributions from K^^rj and Pg lW (P^tt^t] and P^vr^o;) in the final 
state. 


Decay mode 

2P+ - 1 

P/10-^ 

+ 

1 

unknown 

1.09 ±0.16 

? 77 r° 7 r° 

-1 

unknown 

7 ]K+K- 

unknown 

unknown 

rjKpC! 

-1 

unknown 

r^V+TT 

unknown 

0.45 ±0.17 

P 

o 

o 

-1 

unknown 

o 

+ 

1 

-0.288 ±0.036 

unknown 

POP+P- 

-0.136 ±0.088 

unknown 

+ 

1 

o 

o 

unknown 

10.0 ±0.9 

P+P-ttOttO 

unknown 

unknown 

Pg 7 r+ 7 r“ 7 r° 

unknown 

52 ±6 

P° 7 r+ 7 r“ 7 r° 

unknown 

unknown 


unknown 

unknown 

P^TT^TT+TT^Tr^ 

unknown 

42 ±4 

P^ TT^ TT’*' 7 r“ TT’*' 7 r“ 

unknown 

22 ±6 


Treating only vr, P, 7], rj' and u as stable particles, it is easy to see in Table [D that 
almost all two-body CP eigenstate modes have been measured. There are no results for 
the decays > ojrj and ojt]', but these are unlikely to have large branching fractions. 
For three-body decays, in the case that the final state is composed of neutral spin-zero 
particles (Pg l, 7 r°, rj, rj'), a pure CP eigenstate is obtained |5l|. Several such decay modes 
are included in Table [H and there are some others, such as (CP-odd) that are 

worthy of investigation. However, most are unlikely to make a significant contribution to 
the net CP content of P° decays due the cancellation of modes containing Pg^ mesons 
and due to the small branching fractions expected for decays with little phase space. Thus, 
it is unlikely that pure CP eigenstates contribute significantly to Eq. ([9|). 
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In the case of decays to h^h'^h'~ [h^ = Kg^,7r^,ri,ri'] h' = 7r,K), the hnal state 
is in general a mixtnre of CP-even and CP-odd. Specihcally, if the angular momentum 
between the and the h'^h'~ system is L then the CP content is rjho x (—1)'^, where 
rjho is the CP eigenvalue of Hence if both even and odd values of L contribute then 
the hnal state is not a pure CP eigenstate. By conservation of angular momentum, L 
also gives the partial wave in the h'^h'~ system (L = 0 -H- S wave; L = 1 -H- P wave; 
etc.). As discussed above, in the case of 7r’''7r“7r°, the tt+tt" system is dominantly in 
odd partial waves, while for or K^K^K~ both even and odd waves contribute 

approximately equally. For the decays to and it is possible that either 

even or odd partial waves dominate, giving them high net CP content. The dominance of 
even (odd) partial waves is to be expected when the tt+tt" system is mainly isospin zero 
(one), however if the decay is dominated by contributions such as —>• —)■ r^vr^vr^ 

then the net CP content appears hard to predict. In the case, however, of oq dominance, 
a sizable branching fraction for —)■ aQ7r° —)■ (CP-odd) may be expected. 

These modes are certainly worth experimental investigation, as is the as-yet unmea¬ 
sured decay —>■ riK^K~. However, as the branching fractions of the ri7r~^'ir~ 

and r^V+TT” [S3] decays are C(10“^) they cannot make large contributions to the total 
CP content of P° decays. 

As discussed in Sec. O there is a lack of experimental results on hnal states contain¬ 
ing K'l mesons. As the branching fractions for the and K^K^K~ decays are 

relatively high, the residual net CP content between the hnal states containing Kg and 

mesons may be numerically signihcant. In addition to determination of the branching 
fractions, further investigation of the diherences in Dalitz plot distributions between the 
two states with diherent physical kaon states may provide insight into hadronic ehects in 
the multibody decays. 

As the multiplicity of the hnal state increases, the naive expectation is that the net 
CP content is less likely to be signihcantly non-zero, as this would require some form 
of coherence across a multi-dimensional phase space. However, as discussed above the 
hnal state of the 7r+7r“7r+7r“ decay has been shown to be mainly CP-even, so it is 

possible that there are other four- or more-body hnal states with net CP content. The 
hnal states 7r+7r“7r°7r° and appear to be worth investigation; the former is 

known to have a sizable branching fraction [SB]. Similarly to the case discussed above, 
the CP content is given by the angular momentum of the h'^h'~ system, but in this case 
even (odd) partial waves correspond to CP-even (-odd). 

Another four-body decay that is of interest with regard to its CP content is that 
to Pgvr+TT^TT®. The contribution from KgUJ (CP-odd) comprises only around 20 % of the 
total branching fraction, and there may be regions of phase space that are dominated by 
one or the other CP eigenvalue. Although any net CP content would be balanced, to some 
extent, from an opposite ehect in K^7r~^7i'~7r^ decays, there may be a non-negligible 

residual contribution. 

It is likely that a large proportion of the unaccounted-for decay modes are from high- 
multiplicity quasi-flavour-specihc hadronic channels such as K^7r^7i'~^7r~n^, 

iF^7r^7r+7r“7r+7r“, and so on. The net CP-content of such modes is expected to 


be small, as the coherence factors should be low for the reasons discussed above. A great 
deal of additional experimental investigation would be necessary to make a comprehensive 
survey. 

Table [3] contains a selection of modes for which it is not yet possible to evaluate 
the contribution to the reduced width difference. Experimental investigations of these 
channels would be of great interest. 


5 Summary 

There is now a signihcant body of results from studies of various hadronic D meson decays, 
including determinations of net CP content, coherence factors and parameters related 
to strong phase differences. These allow a survey of the contributions to the reduced 
width difference 2 /_d in the system from different types of hadronic decays. If only 

CP-eigenstate and quasi-CP-eigenstates are considered, there appears to be an excess 
of CP-even modes compared to expectation. This demonstrates that the pedagogical 
argument that the width difference arises from such effects is too simplistic. The effects 
of interference between doubly-Cabibbo-suppressed and Cabibbo-favoured amplitudes in 
quasi-flavour-specihc modes are found to be equally important; the small value of yn 
is related to cancellation between the contributions. It would be interesting to see if a 
similar situation occurs in the P°-P° system, where a small reduced width difference is 
expected in the Standard Model m, but enhancements due to new physics effects are 
possible [581 - 160] . Existing measurements of hadronic P° decays are, however, not complete 
enough to allow a data-based study like that presented here; even fewer measurements 
are available in the P° system. 

It has previously been noted that SU(3) symmetry breaking effects can account for 
2 /d ~ 1 % in an “exclusive” approach where D decays are assumed to be dominated by a 
small number of processes [2]. In this approach, the contribution to yn from the SU(3) 
multiplet of decays to hnal states composed of two pseudoscalars is small, but larger 
contributions are expected from pseudoscalar-vector, vector-vector and multibody decays 
as larger SU(3)-breaking effects are induced by phase space considerations. Indeed, it was 
noted that the contribution to yr, from the U-spin doublet of charged K and vr, 

yK. = ^ ^ vr+TT-) + P(P° ^ K+K-) + 2 cos B{D^ ^ K-t:+)) , ( 10 ) 

has almost perfect cancellation; with the beneht of current data it is seen that yxir = 
(0.22 ± 0.30) X 10“^. It is curious to note that of similar relations for modes with an 
additional meson or tt+tt" pair, 

2 /W = 1 ( 2 F+,-l)S(ZlO^/,), ( 11 ) 

i G 7r+7r“7r*^, K~ ^ 

^ (2F+i-l)B(D°-^f,), (12) 

i G 7r+7r“7r+7r”, K'^K~7t'^7t~ , K^7t^7t'^7v~ 
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the former is also consistent with zero within current experimental uncertainties {yK-K-K° = 
(1.7±2.0) X 10“^) while the latter shows a small deviation {yK-Kir+TT- = (—2.8±1.1) x 10“^). 
This may indicate some underlying symmetry of the dynamics of the multibody decays, 
and also appears to support the hypothesis that yo arises largely from SU(3)-breaking 
induced by phase space effects in multibody hadronic D decays. 

Although it is remarkable that existing data on hadronic charm decays allow a survey 
such as that presented in this paper, further experimental investigation would be necessary 
for a full understanding. Improved measurements of the properties of several channels 
are well motivated, including those modes discussed in Sec. S] as well as —>• 

and A'^7r^7r’''7r“. It will also be important to reduce the fraction of unaccounted-for 
decay modes, which are likely to include some high-multiplicity hadronic channels. 
Notably, the BESIII experiment has the potential to improve signihcantly existing mea¬ 
surements, and to add new results, in particular through exploitation of its data sample 
of ^(3770) -7- decays. This sample not only allows challenging decay topologies to 

be reconstructed but also enables measurements that rely on quantum-correlations to be 
performed. Complementary information on hadronic charm decays can also be obtained 
from other experiments such as Belle (II) and LHCb; in addition to branching fraction 
measurements, novel use can be made of charm mixing to determine F+ for quasi-CP 
eigenstates |5], and coherence factors and strong phase differences for quasi-flavour-specihc 
modes | 6 ni 62 ] . 
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